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The paper deals with a planar particle system obeying a generalized exclusion principle (EP) 
' and governed, in the mean field approximation, by a nonlinear Schrodinger equation. We show that 

I the EP involves a mathematically simple and physically transparent mechanism, which allows the 

, ^ ■ genesis of quantum vortices in the system. We obtain in a closed form the shape of the vortices and 

O ' investigate its main physical properties. 

^ ! PACS numbers: 03.65.-w, 03.65.Ge, 05.45.Yv 

In literature experimental evidences has been presented of formation of quantized vortex structures in different 
macroscopic quantum systems. Vortices have been observed both in Bose (Bose-Einstein condensate of alkali atom 
clouds) and in Fermi ('^He-A superfluids, heavy fermion superconductors UPta and Uo.grTho.osBeia) systems [0-^. 
Actually, there exist different theoretical models admitting stationary solutions of vortex type ■ It is apparent 

^ that the vortex type solution is imposed by the non linearities introduced in the particular model adopted. For 
instance, by considering a non relativistic matter system, one can construct different models starting from Non Linear 
Schrodinger Equations (NLSEs) or from the gauged NLSEs in the frame of a Maxwell or Chern-Simons theory. 

Principal goal of the present work is to show that there exists a mathematically simple and physically transparent 
mechanism, imposed by the generalized Exclusion Principle (EP), which allows the genesis of quantum vortices in 
fH , planar non relativistic particle systems. 
Q ' Since 1932, it was clear that the effects due to the statistics and imposed to a system of free fermions by the Pauli 
. exclusion principle can be simulated by a repulsive potential in the coordinate space [T^ ] . In the same way for bosons, 
an attractive potential can simulate the statistical behaviour of the system. After 1940, in different works, particle 
, systems obeying to statistics, which are different from the standard Bosc-Einstcin and Fcrmi-Dirac ones, have been 
> ■ considered [^-|2|. 



, Recently, it has been studied a many body quantum system obeying a generalized Exclusion-Inclusion Principle 
(EIP) p3[-p8||. The peculiar property of this system is that the EIP introduces an attractive or repulsive potential in 
, the coordinate space and can simulate the intermediate (between bosonic and fermionic) behaviour of the system. We 
' consider some examples of real physical systems where EIP can be usefully applied. The Bose-Einstein condensation 
originates from an attraction of statistical nature (Bose-Einstein statistics) among the particles. In several systems 
the Bose-Einstein condensation is studied by means of a cubic NLSE which describes in mean field approximation an 
attractive interaction between two bodies. In place of the cubic and simplest interaction, other interactions can be 
considered as, for instance, the one introduced by EIP to simulate an attraction among the particles. Analogously, in 
the case of superfluids or semiconductors Fermionic systems, the repulsive interaction among the particles originated 
1^ ■ from the Fermi-Dirac statistics can be simulated by the EIP with k < 0. In condensed matter the motion of the 
^ [ couple electron- hole in a semiconductor can be described again by EIP. In fact, while electrons and holes are fermions, 
O ■ together can be considered excited states behaving differently from a fermion or a boson. Moreover, in nuclear physics 
\ the interaction among the fermionic valence nucleons outside the core produces pairs of correlated nucleons that can 
^ ' be approximated as particles with a behaviour intermediate between fermionic and bosonic ones. This nuclear state 

, (quasideuteron state) can be viewed as a particle system that obeys to EIP. 
^ The dynamics of the canonical quantum system obeying the EIP is governed by the following NLSE: 



i^^ = -^^^ + W{p,j)^ + iW{pJ)4' + V4' , (1) 
ot zm 



where the real and imaginary parts of the nonlinearity are given respectively by: 
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The free parameter k is a constant which takes into account the intensity of the exclusion-inclusion statistical effects. 
It is easy to verify that the system described by Eq. ([^) obeys the continuity equation 

| + V.i = 0, (4) 
where p = IV'Pi while the quantum current j is given by: 

j ^-l^{l + np){rVyj-yjVr) ■ (5) 
2 m 

Eq. assures the conservation of the particle number N = J pdPx of the system. 

We discuss now briefly the origin of the model described by Eq. (|^). We start by considering the following nonlinear 
Fokker-Planck equation [E3;E4[ : 



9P „ 



V S 

p(l + Kp) + DV p 



= , (6) 



being 'V S/m the drift velocity and D the diffusion coefficient. In the case k = and D = 0, the above Fokker- 
Planck equation reduces to the well-known continuity equation for p = \ip\'^ and the linear quantum mechanics can 
be obtained if we use the ansatz ■0 = p^/^ exp{i S/h). In Ref. |^ has been considered the case k = 0, D and 
starting from Eq. (H) a new NLSE was obtained. Differently, starting from Eq. (||), after posing k 7^ and D = 0, 
Eq. (|l|) can be obtained. 

The introduction in Eq. (j^) of the factor 1 + Kp originates from the presence of the EIP and allows us to take 
into account many particle quantum effects. In fact, the transition probability from the site x to x' is defined as 
7r(t, X — > x') = r(t, X, x') p{t, x) [1 + k p{t, x')] with r{t, x, x') the transition rate. The transition probability depends 
on the particle population p{t, x) of the starting point a;, and also on the population x') of the arrival point x' . 
For K 7^ the EIP holds and the parameter k, quantifies how much the particle kinetics is affected by the particle 
population of the arrival point. If k > the 7r(i, x x') introduces an inclusion principle. In fact the population 
at the arrival point x' stimulates the transition, and the transition probability increases linearly with p{t, x'). In 
the case k < the 7r(t, x — > x') takes into account the Pauli exclusion principle. If the arrival point x' is empty, 
p{t, x') = 0, and the 7r(t, x — > x') depends only on the population of the starting point. If the arrival site is populated 
< x') < pinax, and the transition is inhibited. The range of values the parameter k can assume is bounded by 
the condition that Tr{t, x x') be real and positive, as r{t, x, x') is. Thus we conclude that k > — l/pmax- 

The form of the nonlinearity W{p,j) in (|l|) is imposed by the continuity equation (^), while the form of W{p,j) is 
imposed by the requirement of the canonicity of the system. The Hamiltonian density of the system ([l]) is given by: 



being 



n = ^\vij\^ + u,,, + vp, (7) 

2 m 

C/b:p=«P^^, (8) 

2 m 

the nonlinear potential introduced by the EIP. In Ref. ||2^ it has been shown that the system (|^) admits 1-D solitons. 

Following the standard procedure of the linear quantum mechanics, we define the quantum velocity v through 
j —vp. Taking into account Eq. and writing ip in terms of the hydrodynamic variables ip = p^^"^ exp(i S/h), we 
have: 

^ VS* 

v = {l + Kp) . (9) 

m 

The quantum velocity can be expressed in terms of the Clebsch potentials S, A, fi through mv = VS' + AV/i. It results 
that the EIP imposes the choice \ = up and p — S for the Clebsch potentials. Now, we define the vorticity uj of the 
system, through a; = V x v and obtain 
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u:^—VpxVS. (10) 
m 

From (p^), it results that Eq. (|l]) describes a vorticose system. In the present contribution, we wiU consider the 
planar, D — 2, static vortex solutions of Eq. (|^) in the case k = — ^ < 0, when EIP is reduced to an exclusion 
principle (EP). By introducing the plane polar coordinates r — \/ x"^ + and 9 = arctan(y/a:) in the plane, we search 
for solutions of Eq. (|l|) in which S — S{6) and p = p{r)- The continuity equation imposes S = HnO and consequently 
we write ^ as: 

V'(r, 0) =p(r)i/2e'"V (11) 

The parameter n must be integer in order to make tp, given by Eq. (ll|), a single value function. In the following, we 
impose / p(Px — N, and therefore p(oo) ~ 0. Let us note that Eq. (|l[) is not Galilei invariant |^8| so that traveling 
solutions can not be obtained by boosting static solutions of Eq. (|^) . 
The quantum velocity v for the vortices (nTl) becomes: 



v = --{l^^p)e, , (12) 
rn r 

where is the unitary vector orthogonal to the vector r = {x, y). After integration of Eq. ( p^ ) on the circle 7^ with 
center at the vortex core and with a radius i? — )■ cxd, we obtain the following relevant property: 

m <j) v-M — 2nh,n, (13) 

which justifies the name of vorticity index n. The vorticity ui ^ lo of the system is given by: 

2nhn . n dp , ^, 

w = 6^ir)-^---f, (14) 

m m r ar 

and taking into account p(oo) — 0, it is easy to verify that the total vorticity depends exclusively on the behaviour of 
the vortex core 

>(fx^2TTn— . (15) 
m 

The definition of the vorticity a; = V x v, imposes the relation: 

uj(fx^ <j> v-d£ , (16) 

which can be easily verified by comparing Eqs (|l^) and (^5|). 

The angular momentum of the vortex, which is a vector orthogonal to the vortex plane, can be calculated as mean 
value of the operator = ~i Tid/dO and assumes the following quantized value: 

L.^nNh , (17) 

which does not depend on the parameter ^. 

We can calculate now the complex non linearity W + iW in Eq. (|l|) for the vortex given by Eq. (^ij) . The non 
linearity becomes a real one and Eq. (|^) reduces to: 

ih^^-—A^-—^pij + V^. 18 
ot 2 m m r'^ 

Eq. ( p^ contains a nonlinearity very close to the one of cubic NLSE and describes a canonical quantum system with 
a Hamiltonian density given by Eq. (^, where the nonlinear potential introduced by the EP assumes the form: 

U..--^^-4p'- (19) 
2 m 



As regard the spatial shape of the vortex p = p{r), we insert Eq. ( pT| ) into Eq. (18), and obtain the following 
second order ordinary differential equation: 
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1 d f dp\ 1/1 dp 



VT 
dr J 2 \p dr 

-^(l-2ep)-^FM=0. (20) 
In the following, we are interested to study the free vortices = 0. After introducing the dimensionless variable 

z = 2n log — , (21) 

where r„ is an arbitrary constant, Eq. ( [20| ) becomes: 

\ ^ 

2^p-l = 0. (22) 



2 d^p n dp 
p dz-^ \p dz 

Now we consider the auxiliary function y{p) 



p dz 

which, taking into account Eq. (E2), obeys the differential equation 



y(p)=(-$^) , (23) 



f + y-l(l-2^p) = 0. (24) 
dp p p 



The function y{p) after integration of (KJ) assumes the form 



2/(p)-- + l-ep, (25) 
P 



being a an integration constant. Combining Eq. ( |23| ) and Eq. (|2^), we arrive to the following first order ordinary 
differential equation for the shape of the vortex: 



dp^ ^ 
dz 



ap + p^l^^p). (26) 




FIG. 1. Plot of the dimensionless quantum velocity 
v' = vmr2/Ji with r2 = yj' C N'/tr as a function of the di- 
mensionless distance from the vortex core r' = r/r2 for the 
vortices with n = 2,4, 6, 8, 10. 



FIG. 2. Plot of the non singular part of the dimensionless 
vorticity ui' = tjmrj/Ti with r2 = ( N/tt as a function of 
the dimensionless distance from the vortex core r' = r/r2 for 
the vortices with n = 2, 4, 6, 8, 10. 
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All the physical solutions of Eq. (^6|), whatever the value of the parameter a € Ft, must be non singular, non 
negative, and normalizable. 

The Hamiltonian H = J T-l<Px of the system can be calculated starting from (Q) and the expression ( p^ for [/^.p! 



H = 



all 

r 



2^ + ::^P^^~^P^ 



(27) 



We must choose a = in order to have a finite value of the Hamiltonian. After integration on the variable 9, we 
obtain the following simple expression for the Hamiltonian 



H = 2TTnh 



jdr 



(28) 



being j = {Tin/mr)p{\ — £_p) the value of the current j — j given by Eq. (^. Finally, also the differential equation 
for the vortex shape assumes a very simple form 



and can be easily integrated, obtaining the vortex profile in an explicit form 

-2 



p{r) = 



Therefore, the wave function of the vortex becomes 

2 



Hr, 0) = 



-\ -1 



exp (i n 6) 



(29) 



(30) 



(31) 



The free parameter r„, which absorbs the integration constant related to Eq. (|29| ), can be calculated from the 
normalization condition 2 tt p{r) r dr — N and assumes the value 



21:1171 , , 8n^ h 
uj^ d^{r) ^ 



2 (n-1) 



(32) 




FIG. 3. Plot of the dimensionless profile § p as a function FIG. 4. 3D representation of the dimensionless profile 

of the dimensionless distance from the vortex core r' = r/r2 for the ground vortex n = 2 in the 2D-dimensionless space: 

with r2 = W (, N/n for the vortices with n = 2, 4, 6, 8, 10. (x' = x/r2, y' = y/r2) with r2 = W ( N/tt. 
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Figure 1 shows the plot of the dimensionless quantum velocity v' = vmr2/fi with r2 = -v/^TV/tt, obtained by 
combining Eqs. ( [l^ ) and (|30|), as a function of the dimensionless distance from the vortex core r' = rjr^ for the 
vortices with n — 2,4,6,8,10. We observe that v is equal to zero for r — r„ where p reaches its maximum value. 
Figure 2 reports the behaviour of the non singular part of the dimensionless vorticity u' = uom r^/Ti, given by (|3^) for 
the same vortices of figure 1. Figure 3 shows the dimensionless profile as a function of r' for the same vortices of 
the previous figures. Finally, a 3D representation of the profile ^ p for the ground vortex n = 2 in the 2D-dimensionless 
space: (x' — x/r2., y' — yli'2) is reported in Figure 4. 



The energy E = H of the n- vortex can be calculated easily by substituting ( |30| ) into (28) and performing the 
integration 

E-\n\^ ; k^^- . (33) 

The energy of the system results to be quantized and the energy spectrum lower bounded. 

We recall that the family of the vortices (30) is related to the vortices Pjp{r) of ref. through p{r) cx r'^Pj^ir). 



We remark that the vortices Pjp [r) are obtained as self-dual static solutions of a Chern-Simons model and correspond 
to the same energy state with energy equal to zero. On the contrary, for the vortex family (^), we have that any 
vortex corresponds to a different state of the system whose energy is given by Eq. (|3^) . 

We conclude by noting that very recently |^ , it has been considered a generalization of the present model (describ- 
ing neutral particles which obey the EIP) in the cases of non relativistic charged particles. In this modified model the 
matter field obeying the EIP is minimally coupled to a gauge field whose dynamics are described within the frame of 
the Chern-Simons picture. The model is a canonical one and admits self-dual static non topological vortex solutions 
with zero energy and linear momentum. The expressions of the main physical quantities associated to these solutions 
are obtained. The electric charge and the angular momentum are derived analytically, while the shape, together with 
the electric and magnetic fields of the vortex, are obtained numerically. This model can be considered as a continuous 
deformation of the Jackiw and Pi one ||l^ , performed by the parameter n which takes into account the EIP. 
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